Coarse grained and fine dynamics in trapped ion 

Raman schemes 



B. Militello 1 * P. Aniello 2 ! and A. Messina 1 1 

1 Istituto Nazionale di Fisica della Materia, Unita di Palermo, MIUR 
and Dipartimento di Scienze Fisiche ed Astronomiche dell'Universita di Palermo, 
Via Archirafi, 36 - 90123 Palermo, Italy 
Istituto Nazionale di Fisica Nucleare, Sezione di Napoli, 
and Dipartimento di Scienze Fisiche dell'Universita di Napoli "Federico II" , 
Complesso Universitario di Monte S. Angelo, Via Cintia - 80126 Napoli, Italy 

February 1, 2008 



Abstract 

A novel result concerning Raman coupling schemes implemented using 
trapped ions is obtained. By means of an operator perturbative approach, 
it is shown that the complete time evolution of these systems can be ex- 
pressed, with a high degree of accuracy, as the product of two unitary 
evolutions. The first one describes the time evolution related to an effec- 
tive coarse grained dynamics. The second is a suitable correction restoring 
the fine dynamics suppressed by the coarse graining performed to adia- 
batically eliminate the nonresonantly coupled atomic level. 
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1 Introduction 



Trapped ions provide an effective platform for observing interesting aspects 
of quantum mechanics and for realizing useful applications in the context of 
quantum computation mESEJ. 

In these traps, a time-dependent quadrupolar electromagnetic field is re- 
sponsible for a charged particle motion which may be kinematically assimilated 
to the motion of a massive spot subjected to a quadratic potential. Such a 
circumstance provides the possibility of describing the center of mass of an ion 
confined into a rf Paul trap as a quantum harmonic oscillator . In addi- 
tion the ion possesses atomic degrees of freedom related to the electronic motion 
around the nucleus H 03 El 

Acting upon the system via laser fields, it is possible to induce vibronic 
transitions described by Jaynes-Cummings-like Hamiltonians |S| . Such in- 
teractions are characterized by nonlinearities governed by the so called Lamb- 
Dicke parameter which, in a spherically symmetric trap, is nothing but the ra- 
tio between the width of the ion vibrational ground wave function and the laser 
wavelength. Controlling the Lamb-Dicke parameter leaving unchanged the laser 
frequency would then provide the possibility of implementing a wider variety of 
Hamiltonian models. Unfortunately, these two parameters are strictly related, 
the Lamb-Dicke parameter being proportional to the inverse of the wavelenght, 
hence proportional to the laser frequency. Therefore, they cannot be indepen- 
dently adjusted. Nevertheless, such a possibility can be achieved exploiting a 
fundamental property of Raman coupling schemes. In these couplings a three- 
level system is subjected to a far detuned A-scheme. Under some assumptions 
and within approximations concerning the time scale used to observe the sys- 
tem, the dynamics may be well enough described via an effective Hamiltonian 
thinkable as a Jaynes-Cummings-like Hamiltonian related to an effective laser 
field, with frequency and wave- vector given by the differences between the corre- 
sponding parameters of the two real A-scheme lasers. Then, changing the angle 
between the two real laser propagation directions leads to the possibility of ob- 
taining effective lasers such that the product of their wavelenght and frequency 
is not the velocity of light. The price to pay is the complete ignorance about 
the detailed system dynamics at a fine time scale. 

In this paper we will try to overcome such a limit by means of a new ap- 
proach. More precisely, we will analyze the dynamics of a three-level trapped 
ion subjected to a Raman A coupling scheme using a very convenient perturba- 
tive decomposition of the evolution operator of the system. We will show that 
the Raman scheme time evolution can be factorized, at the second perturba- 
tive order (hence with a high degree of accuracy), into two unitary evolutions. 
The first one can be interpreted as the effective time evolution which may be 
obtained adiabatically eliminating one of the three atomic levels involved in the 
coupling scheme Such a dynamics concerns the coarse grained variables and 
is in accordance with already well known results [2] • The second unitary evolu- 
tion introduces the correction necessary to take into account the fine deviation 
from the coarse grained time evolution. 
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2 The physical system 



The physical system on which we focus is a three-level harmonically trapped 
ion subjected to a Raman coupling scheme involving atomic transitions. The 
relevant Schrodinger picture Hamiltonian of the Raman A-scheme is then given 

by 



1=1,2,3 a=x,y,z 

+ [frgiae-^r-^ata + h.c] + [hg 23 e-<^-^ f ) a 23 + h.c] ,(1) 

where au = \k) (l\ (with k,l = 1,2,3), {\l)} being the considered three atomic 
levels and {hioi} the corresponding energies; a a (a = x,y, z) is the annihilation 
operator related to the center of mass harmonic motion along the direction a 
(we will denote the associated Fock basis by {ip"})- For the sake of simplicity 
(but without loss of generality), we have assumed to deal with a 3D degenerate 
parabolic trap with single frequency v. The two laser fields responsible for the 
coupling terms are characterized by complex strengths (proportional to the laser 
ampliutude and to the atomic dipole operator, and including the laser phases), 
wave vectors and frequencies g 13 , fc 13 , w 13 and g 23 , ^23, ^iz respectively. 

The level |3) is assumed to be dipole-coupled to both the levels |1) and |2) 
via far detuned lasers. Precisely, the two laser frequencies are chosen in such a 
way that 

A = uj 3 - u>i - lo 13 = uj 3 - uj 2 - UJ23, (2) 
where the detuning A satisfies the condition 

|A| » IsiaMsasI,!'. (3) 

The analysis of such an Hamiltonian model has been already carried out, 
for instance in ref. providing the adiabatic elimination of the nonresonantly 
coupled atomic level |3) following the path pointed out in refs. ^3 and 
including the motional degrees of freedom. Indeed, due to the large detuning, 
the transitions for instance from the level 1 1) to the level |3) are very fast and 
immediately followed by decays on the atomic level |2). Therefore, considering 
only coarse grained observables, meaning that the system is observed at a "rough 
enough time scale" , effectively eliminates the far detuned level; namely, at such 
a time scale, the only observables and hence meaningful dynamical behaviors, 
involve levels |1) and \2) as a result of time averaging second order processes 
having |3) as an intermediate virtual level. This procedure hence suppresses the 
fine dynamics, that is it sacrifices any information concerning the fast dynamics 
the third level is involved in. 

In the following we present a perturbative approach to the solution of the 
dynamical problem related to Ha overcoming the limit of the coarse graining, 
making it possible to study also the fine dynamics discarded by the adiabatic 



3 



elimination. The first step consists in passing to a rotating frame, meaning that 
the time-dependent Hamiltonian H\ is canonically transformed via the operator 



R(t) = e - tAt , (4) 

where 

A = — A) (<Tn + 022 + O33) - Wi 3 (7n — W23CT22 

= w\au + W2CT22 + (<^3 - A) ct 33 , (5) 
into the following time-independent rotating frame Hamiltonian: 

h ■= R(tV (-HaW - 4) = frAH = hA (h + H b + flj) , (6) 

where H is a dimensionless Hamiltonian which is the sum of the three hermitian 
operators Hq, Hb and i/j defined as 

= 0-33, 

_ A 2^a=2;, a ,2 a a a ai (7) 



Considering the assumption given by the inequality both £fe and /ft 
may be thought of as perturbations with respect to Hq. In fact, introducing the 
dimensionless perturbative parameter 

A:= X' 3 = m ax{^, |gi 3 |, I523I}, (8) 

both Hb and fit are first order perturbations in A: 

H = H(X) =ffo + Ax ^ ala a + A ^ [x i3 e^' 3 '^ + /i.e.] , (9) 

ot—x,y,z j — 1,2 

where k = vjg < 1, Xj^ = gj^/g, < 1, and we notice that, due to 

condition @, Acl. It is worth noting that the circumstance that Hb is treated 
as a perturbation leads to the eccentric situation of an unperturbed Hamiltonian, 
Hq, wherein the bosonic degrees of freedom are absent. Nevertheless, as we shall 
see, such a mathematical artifice reveals fruitful in order to succeed in factorizing 
the coarse grained dynamics and its fine correction. 

Our solving procedure relies on a suitable canonical transformation e %z ^ of 
the rotating frame Hamiltonian such that 

e i*W e -&m = h A e l2{X) H{\) e~^ (A) = HA (h + C*(A)) , (10) 

where C(X),Z(X) depend analytically on the perturbative parameter A and 
C(A) is a constant of motion with respect to the unperturbed dynamics, i.e. 
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[Ho,C(X)} = 0. This transformation allows to give a very convenient decompo- 
sition of the evolution operator associated with the rotating frame Hamiltonian, 
namely 

(id 

At this point, truncating the power expansions 

(7(A) = XC 1 + X 2 C 2 + ■■■ + X n C n + ■ • • , Z{X) = XZ l + X 2 Z 2 + ■■■ + X n Z n + ■■■ 

at a given perturbative order, one obtains by formula (fill) useful expressions of 
the evolution operator. This procedure has been developed in a general setting 
in refs. ^JEj- In the next section, we want to recall briefly the mathematical 
background and to show how the operators C\, C%, . . ., Zi, Z%, . . . can be com- 
puted by a suitable iterative process. Then, we will give the explicit solutions 
for our case up to the second perturbative order. 

3 Perturbative analysis of the rotating frame 
Hamiltonian 

Let H u , Hp be hermitian operators and assume that H u has a purely discrete 
spectrum. Denote by E < E\ < E 2 < ■ ■ ■ the (possibly degenerate) eigenvalues 
of H u and by Pq, Pi, P 2 , ■ ■ ■ the associated eigenprojectors. Now, consider the 
operator H(X) = H u + XH p , X G C, which is hermitian if A is real. It is pos- 
sible to show that, under certain conditions J3Jj there exist positive constants 
ro,r\,r 2 , . . . and a simply connected neighbourhood X of zero in C such that 
the following contour integral on the complex plane 

An (A) = ^-<f dE (H{X) - E) A e X, 

llx J\E-E m \=r,} ' 

defines a projection (P m (X) 2 = P m (A)), which is an orthogonal projection for 
real A, with P m (0) = P m , and IsAh P m (X) is an analytic operator- valued 
function. Moreover, the range of P m (A) is an invariant subspace for H(X), hence 

H(X) P m (X) = P m (X) H(X) P m (A), (12) 

and there exists an analytic family U(X) of invertible operators such that 

P m = L>(A) J P m (A)L>(A)- 1 , U(0)=ld, (13) 

and U(X) = e l2(x \ X E 1, with Z(X*) = Z(A) t (hence, for real A, Z(X) is 
hermitian and U(X) is unitary), where IbAh Z(X) is analytic. One can show 
easily that the function A ^ U(X) is not defined uniquely by condition l|13fl 
even in the simplest case when H u has a nondegenerate spectrum. Anyway, the 
nonunicity in the definition of U(X) is not relevant if one is only interested in 
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obtaining an expression of the evolution operator associated with H(X) of the 
general form We will see soon that there is a natural condition which fixes 
a unique solution for U (A) . 

Now, let us define the operator 



K(X) := U(X)H(X)U(X)-\ 



(14) 



which, for real A, is unitarily equivalent to H(X). Using relations l(T2|) and (|T3|) . 
we find 

K(X) P m = U(X) H(X) P m (X) U(X)- 1 = U(X) P m (A) H(X) P m (X) U(X)- 1 

and hence: K(X) P m = P m K(X) P m . It follows that H u , K(X) = and then 
we obtain the following important decomposition formula: 



U(X) H(X) UiXy 1 = H u + (7(A), 



(15) 



where 



C(X),H n = 0, i.e. C(A) is a constant of the motion with respect to 

the time evolution generated by H u . At this point, we can obtain perturbative 
expressions of the unknown operators C'(A), U(X) by means of a recursive alge- 
braic procedure. 

Indeed, since the functions A <—> C(X) and A <—> Z{X) are analytic in X and 
(7(0) = Z(0) = 0, we can write: 



(7(A) = A" C n , Z(X) = J2 xn Zn AG J. 



(16) 



n=l 



In order to determine the operators {C n } and {Z n }, we substitute the exponen- 
tial form e lZ ( x ^ of U(X) in formula (|15|l thus getting 

H(X) + jr ^ ad| (A) H(X) = H u + (7(A). 



n=l 



where we recall that ad^ A ^-ff(A) := [Z{X),H{X)\. 

Next, inserting the power expansions (|16(l in this equation, in correspondence 
to the various perturbative orders, we obtain the following set of conditions: 



C 2 



Zi , H u 



l r 
2 



Zi, 



Zi,H u 



Zi,H u 



Zi,H p 



Ci,H u 
C*2 . H u 



where we have taken into account also the additional constraint [(7(A), H u ] = 0. 
This infinite set of equations can be solved recursively. The first equation, to- 
gether with the first constraint, determines Z\ up to an operator commuting 
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with H u and C\ uniquely and so on. It is convenient to eliminate the arbitrari- 
ness in the determination of the operators {Z n } choosing the minimal solution 
characterized by the additional condition ^2 m P m Z n P m = 0, n = 1, 2, . . . . 
In our particular case, we have the following identifications: 



H u = H , 

Hp = Hb + Hi ■ 



(17) 



Notice that the two (infinitely degenerate) eigenspaces of the unperturbed Hamil- 
tonian H are associated with the eigenprojectors 



{P m } m =g,e = {Pg = lfi ® (<Tll + CT22) , P e = Ifi ® (£33)}, 



(18) 



where 

E (NO^»J)®"-®(|VO(^nJ) (19) 

is the identity in the vibrational Hilbert space. Accordingly for the operators 
{Ci, Z\, C2, Z 2 , . . .} forming the minimal solution, we get at the first perturba- 
tive order the following expressions: 



ACl = E ro =e, S ^(^ S +^l)4 

Similarly, at the second order, we have: 



^2 = Em=e, g Pm{i[Zl,H B +Hl 

>?Z 2 = iY,^ k Pi{i[Zi,H B + Hi 



Z^Z^H^Pk. 



Eventually, performing explicit calculations, we find that 

ACi = Hb, 



(20) 



(21) 



(22) 



x2A |gl3| 2 
C2 



\g23W , |gi 3 | 2 + |g 23 | 2 . 

Cll 1^5— ^22 H 7^5 CT33 



/ 513532 „—ik 13 -r e ik23-rfi.^ _|_ 



A 2 



A 2 
/i.e.) , 



V A 2 

where we have set 33 j = 5^3, arid 

AZ a = 1 ^e-^ p a 13 - h.c) + 1 (^e"* 3 ^ - h.c) , 



(23) 



(24) 



A 2 Z 2 = I ^^-^"13(713 + ^X 3 ia-3i^ + (^-X 2 3<723 + ^-£32032) | , (25) 
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where: 



\ X 3j := i[e^, Z a=x , y , z 4«J = Xj 3 , 

with j = 1, 2. 

The interpretation of this result leads to a very interesting fact. Indeed, it 
turns out that once the unitary transformation e lZ ^ has been applied to the 
rotating frame Hamiltonian S) (recall eq. (|10p). the time evolution of the system 
is described, at the second order in the parameter A, by the Hamiltonian 



HA (ff + Ad + X 2 C 2 j = £12 + Sj 3 , (26) 
where, in order to display a more transparent formula, we set 

f)l2 := Hv 22 (®a®a) ® (^11 + ^22) + ^WliTll + hw 2 U 22 

a=x,y,z 

+ \hg 12 e- l ^a 12 + h.c. 



(27) 



fls := V (oj, o a ) ® <r 33 + fi(A + w 3 )a 33 , (28) 



with: 



|3j3| 2 . . „ . I.9i3 1 2 + 1^23 1 2 513332 r r r 

—1—) J = I? 2 , cj 3 = , 312 = — - — , fcl2 = fcl3 - fe 23 - 



Thus, the transformed Hamiltonian is the sum of two decoupled Hamiltonians 
f)i2 and ,^3, [#12,^3] = 0, 'living' respectively in the ranges of the orthogonal 
projectors P g and P e . This is a consequence of the fact that [P m ,C n ] = 0, 
to = g, e , n = 1, 2, . . . . It is worth noting that the Hamiltonian ,f)i2 can 
be regarded as the rotating frame Hamiltonian of a trapped two-level ion in 
interaction with a laser held characterized by the following parameters: 

r CJ12 = CJl 3 — UJ 32 = L0 2 — (J I, , , 

lfci2= k 13 -k 23 . 

This effective coupling can be compared with the result found performing the 
adiabatic elimination of the level |3) (see ref. 9 ). We will come back to this 
point in the next section. 



4 Dynamics of the Raman scheme 

The question of what the complete dynamics of the system is now arises. First, 
it will be convenient to adopt the following notation. Given a couple of functions 
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/ and h of the perurbative parameter A, if /(A) = h(X) + 0(A 3 ), we will write 
simply: 

/(A) ~ h(X). 

Next, let us denote by T\ the evolution operator associated with the Raman 
scheme: 

%K {jt fA ) (<) = ^ A (t} ^ A(<) ' ^ a(0) = Id ' (30) 

Expressing T\ in terms of the evolution operator associated with the rotating 
frame Hamiltonian yields: 

T A (i) =R(t) exp(~htj . (31) 

Now, according to what we have shown in the previous section, we have: 

f(t) := exp(~ht 

= exp(-zA H(X) t) 

£ e -j(AZ! + A 2 Z 2 ) e -iA(H +\C 1 +\ 2 C 2 )t e i(XZ!+\ 2 Z 2 ) /g 2 ) 

where we have truncated the power expansions of Z(X) and C(X) at the second 
order in A. Formula 1|32|) provides an approximate expression of the evolution 
operator in the remarkable form of a one-parameter group of unitary transfor- 
mations. Nevertheless, in order to achieve an approximate expression allowing 
a direct comparison with the coarse grained dynamics, we still need to per- 
form some manipulation. To this aim, observe that, since the commutators 
[A 2 C 2 , AZi + A 2 Z 2 ] and [X 2 C 2 , H + ACk] = [A 2 C* 2 , XC{\ are of the third order 
in A, maintaining our degree of approximation we can write 

e -i(XZ 1 + X 2 Z 2 ) e -iA(H + XC 1 + \ 2 C 2 )t £ e -i(AZi + A 2 Z 2 ) e - l AA 2 C 2 t e -iA(H +ACi)t 

£ e --iAA 2 C 2 t e -i(AZ 1 + A 2 Z 2 ) e -iA(fl +AC' 1 )i 

Therefore, we can manipulate the second order expression i|32|) of T(t) as follows: 

f(t) « e - iAA2 C , 2t e -j(Aii+A 2 Z 2 ) e ~jA(H +AC'i)t e i(Ai 1 +A 2 Z 2 ) 
. e -jAA 2 C 2 t e -iA(H +AC'i)t 
x e iA{Ho+\C 1 )t e ~i{\Z 1 +\ 2 Z 2 ) e -iA{H +\C 1 )t e i{\Z 1 +X 2 Z 2 ) 

Finally, we find a remarkable decomposition of T: 

f(t) « f e (t)ff(t), (33) 
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where we have set 

f c (t) := exp(-iA (h + \& + A 2 C* 2 ) tj , (34) 

f f (t) := exp(-i (AZi(t) + A 2 Z 2 (t))) exp(z (aZi + A 2 Z 2 )) , (35) 

with Z k (t) = e iA(H +ACx)t £ fe e - 4 A(ff +AC 1 )t ; fc = 1; 2. it is worth emphasizing 
that, by a completely analogous procedure 1 , we get also 

f (t) « f/(t) T (f), T/W := f f (-i) f - (36) 

Notice that, due to the specific dependence of Tf on t, one has that T f '(t) ^ Tf (t). 

In the light of formula ilMMI) (or 1)36(1). the time evolution in the rotating 
frame, given by T, may be thought of as a process consisting of two fundamen- 
tal components. One of these is an effective time evolution, described by T e , 
in which the levels |1) , |2) are decoupled from the level |3). The other compo- 
nent, the one Tf (or T f ') is responsible for, is a correction to T and involves fast 
transitions (the operators Z\(t) and Z%{£) oscillate at the detuning frequency 
A) from and to the third atomic level. Considering the complete time evo- 
lution H31f) . observe that the unitary evolution described by RT e corresponds 
to the effective dynamics obtained in [5] restricting the analysis to the coarse 
grained observables. In fact, RT C is the evolution operator associated with the 
time-dependent effective Hamiltonian 

H e = HW+H&, (37) 

where 

H^ 12 \t) := TlV ^ (4«a) ® (<5ll + ^22) + h((Jl + #1)0-11 + h(ld 2 + W 2 )<T22 

a— x,y,z 

+ \hg 12 e~ l ^~^d 12 + /i.e.] , (38) 
:= hv J2 (SaO® *ss + h(u> 3 + #3)0-33. (39) 

We remark that we have deduced this result analytically; no adiabatic approx- 
imation has been performed. We also stress that a relevant difference between 
T and Tf consists in the kind of time dependence. Indeed, on one hand, the 
unitary evolution T e forms a one-parameter group, hence it is expressible as the 
exponential of a generator multiplied by t. It follows that a truncated power 
expansion of the exponential retains its validity only on a finite time span. On 



1 The calculation may be carried on directly, i.e. step by step as in the previous case but 
changing the reordering of the exponentials in 1321 and subsequent formulae, or from 1331 
exploiting the fact that f(t) = f(-t) t . 
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the the other hand, Tf can be expressed as the exponential of an operator whose 
time dependence involves only sinusoidal factors. In fact, we have: 

T f « e - i ( A (^i( t )"^i)+ A2 (^(«)-^2)) e ±[\z 1 (t)+\ 2 z 2 (t),\z 1 +\ 2 z 2 ] 

£ e -i{X{Z x {t)-Z 1 )+\ 2 {Z 2 {t)-Z 2 )) e ±\ 2 [Z 1 (t),Z 1 ] 
_ p -i(A(Zi(t)-ii)+A 2 (^2(t)-i2)) 

where \Z\(t),Z\\ — has been used. It then follows that the truncated expan- 
sion 

T « 1 - zA (^(t) - Zf) - i\ 2 (Z 2 {t) - Z 2 ) - h 2 (Zi(t) - if)' (40) 

is legitimated independently on time. 

Summarizing, we have shown that the standard adiabatic elimination tech- 
nique provides an effective dynamics, described by RT e , that differs from the 
complete second order dynamics of the Raman scheme for the presence of an- 
other unitary evolution which can be cast in the form of the exponential of a 
rapidly oscillating operator function of time. Therefore, the factorization into 
a coarse grained and a fine dynamics given by eq. ll-i-il) makes the correction to 
the adiabatic approximation solution very readable and easy to be calculated, 
in view of the expression (|40|l . It is worth noting that the corrections due to 
Tf are small in amplitude, since the operators XZ±, \Zi(t) and A 2 Z 2 , A 2 Z 2 (t) 
are respectively of the first and second order in the perturbative parameter. 
Moreover, they provide terms oscillating at the detuning frequency. Hence, as 
expected, the fine dynamics is small and fast. As a conclusive remark, we wish 
to emphasize that this micro-fast behavior brought to the light with the help of 
the method applied here should become of practical interest in connection with 
time resolution improvements of experiments. 
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